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Abstract 

We propose using an array of potential wells as an interferometer, in which the beam splitters are 
provided by tunneling during an appropriate time through the barrier between wells. This arrange¬ 
ment allows demonstration of generalized Hong-Ou-Mandel effects with multiple particles traversing 
one or several beam splitters. Other interferometer effects can occur, including a violation of the 
Bell-Clauser-Horne-Shimony-Holt form of the Bell inequality. With interactions, one sees various 
effects including so-called fermionization, collective tunneling, and self-trapping. 
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1 Introduction 

The Hong-Ou-Mandel (HOM) effect [1, 2] is perhaps the simplest example of boson-boson interference, 
in which two particles approaching in two different inputs of a beam splitter always exit both in the same 
output, because destructive interference removes the possibility of exits in separate outputs. Experimental 
verifications have been carried out with photons [1, 3]. Other cases with two photons from each side [4] 
as well as two from one side and one from the other have also been considered theoretically [5], and 
experimentally seen [6]. Theoretical generalizations using cold atoms have recently been considered, 
including one in which Fock states of arbitrary numbers of particles impinge on a beam splitter [7, 
8], and another using pair-correlated atoms produced via a collision of two Bose-Einstein condensates 
[9]. Recent experiments using particles rather than photons have involved electrons [10], ®^Rb atoms 
trapped in optical tweezers [11], and helium atoms [12] where the beam splitter was an Bragg scattering 
optical grating. Each of the papers using matter waves rather than photons emphasizes that experiments 
developing coherent indistinguishable pairs of particles may be relevant in other areas such as quantum 
computing and information processing [13], highly sensitive force detection [14], quantum simulations 
[15], testing Bell inequalities with material observables [16], etc . Here we want to consider HOM-like 
interference with the use of Fock states of cold atoms undergoing simple tunneling in optical double (or 
multiple) potential-wells in a method analogous to that experimentally used in Ref. [11]. 

The essential element in the HOM effect, or in any interferometer, is a beam splitter. In the case of 
cold gases, various forms of beam splitter have been designed, including, e.g., Bragg scatterers [12, 17] 
and double potential-well devices. In the latter case, the usual format [18] involves guiding a matter wave 
through a potential well having a spatially or temporally growing central repulsive peak that divides the 
beam into two parts, which, for example, can be recombined further along to show interference effects 
analogous to a two-slit interferometer. In one case the authors [19] proposed using a pair of side-by-side 
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wave guides having a narrow region where tunneling took place; by adjusting the wave packet width (or 
momentum) and the barrier width the packet splits into two equal parts, thus using the time of tunneling 
as the main element to form a 50-50 beam splitter. Recently Daley et al [20] suggested using double-well 
tunneling as a beam splitter to measure the order-2 Renyi entropy of an entangled state. Compagno et 
al [21] considered boson and fermion walkers on a one-dimensional lattice with an optical impurity with 
transmission and reflection coefficients adjusted to act as a beam splitter. Here we will further examine 
the possible use of the time of tunneling as a beam splitter, a method already adopted by Ref. [11] in 
their experiment. 

Most interference experiments with cold atoms have involved observing a periodicity in the particle 
density as a function of position or angle. An alternative would be counting the number of particles in 
the detectors placed after a beam splitter. In the case of the photon HOM experiment a minimum in 
the coincidence between the two detectors was observed. However, there has been considerable recent 
progress in the actual observation of individual atoms via the “quantum gas microscope” [22, 23[. In 
this case the “approach is to assemble quantum information systems with full control over all degrees of 
freedom, atom by atom, ion by ion” [22]. We will assume such an approach is literally possible so that the 
number of atoms in each well is known at the beginning of the experiment (input to the beam splitter) 
and at the end in the detectors after the beam splitters. As in Ref. [11] we assnme that the beam splitter 
involves turning on tunneling between two or more wells for a set time. In this experiment two atoms 
were trapped in their ground states in wells with optical tweezers and observed via fluorescence. The 
distance and barrier between the wells could be varied to enhance or cut-off tunneling. In Ref [12] Bragg 
scattering was used for mirrors and a beam splitter with observation of pairs or singles of atoms falling 
on a mnltipixel detector. We will look at variations of the HOM effect but also see that these remarkable 
experimental methods may yield other interesting possibilities where more elaborate interferometry is 
possible. 


2 Double-well beam splitter 


Fig. 1(a) shows a beam splitter with two particle sources and two detectors. In the standard beam 
splitter with the annihilation operator for source A being a and that for B is (3 then the annihilation 
operators for detectors 1 and 2 are, respectively 


di 

d2 


-^{a + il3) 

T(.a + m 


( 1 ) 


We want to show that the the double well of Fig. 1(b) is equivalent to the device of 1(a) if the 
tunneling occurs for a set time period. Let a be the destruction operator for a particle in well A and b 
for one in well B. Then the two-level Hamiltonian we consider is taken as 

H/h = i?o(a^a + — \{a^b -\- b^a) -I- —{a^a^aa + b^b^bb) (2) 

where we have assnmed the wells are at equal depth with identical particles in each. The tunneling 
parameter is A with interactions of strength W among particles in the same well. In the double-well case 
with no interactions the annhilation operators evolve in time according to 


a{t) = a{0) cos \t + ib{0) sin \t 
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Figure 1: a) Na,Nb bosons proceed from the sources A and B to a beam splitter, followed by two 
detectors 1 and 2, where mi and m 2 particles are detected, b) A double-well potential with tunneling 
strength A between cells A and B. 


b{t) = ia{0) sin Xt + b{0) cos Xt (3) 

At time t = tt/AX we have 

“(ix) ^ (a(0) + *^(0)) = 

^ (*«(0) + ^(0)) = ^2 (4) 

which is then equivalent to the 50-50 beam splitter. The simple two-particle HOM effect, with a single 
particle initially in each well, corresponds to the following state vector: 

(ix) (ix) = X = T + bioybio)^) |0) (5) 

that is, ending up with a superposition of a pair of particles in each well as expected in the simplest HOM 
effect. The HOM effect arises because of the destructive interference resulting from exchange of particles 
in the final state. 

For the general case of arbitrary N we can solve for the coefficients in the expansion of the wave 
function |■i/'(^)) of ^ Fock states |n,iV — n), having n particles in well A and N — n in well B. If we 
take 

li’it)) ='^Cnit)\n,N-n) (6) 

71 

then the differential equation for c„(t) can be shown to be 

i — {n,N — n\H/h\n',N — n') Cn'{t) 
at ^ 

n' 

= -X |^\/n(A^ -n + l)c„_i(t) + ^/{n+ 1)(A - n)c„+i(t) 

+ ^ [n^ - n + {N - nf - {N - n)\ Cn{t) (7) 

where we have omitted the term in Aq that leads to a trivial phase factor. 
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3 HOM calculations 


We can use Eqs. (7) to compute HOM-type interferences. For the simplest case of Na = Nb = 1, the 
equations are 

ico = jco — V2ci 

ici = —V2co — V2 c2 ( 8 ) 

ic2 = 7 C 2 — V^ci (9) 

where now we have expressed the time in units of fi/A and defined 7 = WjX. The solutions of these 

equations for the probabilities |c„(t)p are shown in Fig. 2 

k„(t)p 


t 

Figure 2: Double-well HOM result with N = 2, Na = Nb = 1, showing the time dependence of the 
probabilities of particle occupation of the wells with no interaction, 7 = 0 . The probabilities are shown 
with |ci(t)p as a solid line, and |co(t)P = |c 2 (t)P, dashed. Initially we assume the wells each contain one 
particle: ci(0) = 1, and co(0) = 02 ( 0 ) = 0. After time t = 7r/4 (in units of h/X) the probabilities become 
jcip = 0 , and |coP = |c 2 p = 1 / 2 . 

For 7 = 0 the analytic solutions are 

'I 

co(t) = C 2 (t) =-p sin2t 
v 2 

ci(t) = cos 2 t ( 10 ) 

giving the wave function at time t = tt jA as 

|V^(t = 7r/4)) = ^(|0,2) + |2,0)) (11) 

Thus the probabilities for the HOM effect, with 

= |c„( 7 r/ 4 )p (12) 

are the standard results: po = P 2 = 1/2; pi = 0. 
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This “HOM time” is just as identified in the experiment of Ref. [11], but in our case, the HOM effect 
can be generalized to more particles in each source as was done in Refs. [7, 8 ]. In Fig. 3 we show the 
cases of iV = 4 and 8 , where again we consider the occupations of the wells at time t = 7 r /4 with equal 
occupation in each well initially. If the initial number of particles is not the same in each well then the 


Pn Pn 



Figure 3: Double-well HOM result with = 4 and 8 . (Left) Probabilities of particle occupation of 
the well A at time time t = 7 r/ 4 , for initial conditions Na = Nb = 2; the interaction parameter 7 = 0 . 
(Right) Similarly for initial conditions Na = Nb = 4. 

HOM preference for even occupation can’t be maintained so purely as seen in Fig. 4. 

Pn 



Figure 4: Double-well HOM result for N = 9, with Na = 4, Nb = 5 initially, showing the probabilities 
Pn of particle occupation of the well A at time time t = 7 r/ 4 . 

Until now we have reproduced the results of Ref. [7] exactly. One aspect of the present approach to 
the HOM analysis is that we can include the effects of interactions. In the N = 2 case the probabilities, 
gotten from a numerical solution of Eqs. ( 8 ), show that po and p 2 decrease as 7 increases and pi increases 
until all probabilities are equal at 7 ~ 2.5. By 7 = 6 , po ^nd p 2 are very small. This might be expected 
for a repulsive potential, but the result, (and the Eqs. (7) for any time) are invariant under sign reversal 
of 7 . The results are the same as if the particles had been fermions. Indeed such so-called fermionization 
of the system has been seen previously [21, 24]. We can see what is happening by looking at the full time 
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k„(t)p 



Figure 5: Double-well HOM probabilities |c„(t)p as a function of time for = 2 (initially Na = Nb = 1) 
with 7 = 6 ; |co(t)P, |c 2 (t)P shown as dashed lines; |ci(t)p as solid line. Compare with Fig. 2, for which 
7 = 0. As 7 is increased beyond the value here, the amplitudes of the oscillations diminish further. 


dependence in Fig. 5. The average probability for staying in the original {1,1} state is enhanced and 
those for the pair states {2,0} and {0,2} are diminished. This result illustrates coordinated tunneling as 
seen in other work [24, 25, 26]. In Refs. [24] and [26] two particles were started out in the same well, and 
the result was the pair states became favored and the {1,1} state diminished in probability. In our case 
any coordinated tunneling would seem to involve the pairs crossing together, but in opposite directions, 
so that we get the opposite result. In Ref. [21] bunching was seen for two non-interacting bosons and 
antibunching for a pair of fermions or strongly repulsive bosons. 

In Fig. 6 we show the effect of increasing 7 = WjX on the initially equal-sided N = 8 case. At small 
7 = 0.3 the HOM effect becomes less “pure” with odd states {1,7},{3,5}, etc now possible; by 7 = 0.5 
the {4,4} state is enhanced; and at large 7=1 the {4,4} state has become highly favored. The effect 
seems related to “self-trapping” in which the sign of an initial imbalance in the number of particles in 
the two wells becomes fixed when the interaction parameter exceeds a critical value [26, 27, 28, 29, 30], 
with reduced amplitude oscillations away from the initial value of the imbalance. We discuss some results 
of an analysis of self-trapping in the Appendix and how onr results might be similar. For fixed N = 8 
the critical interaction strength for self-trapping is 7 = 0.5; the graphs then show results below, at, and 
above this critical value. This effect is again independent of the sign of 7 . 

4 Generalizations 

In addition to simulating a single beam splitter, particle tunneling between wells can behave as an 
interferometer. Ref. [31] coupled two wave guides together with tunneling at two points, corresponding 
to beam splitters, to make a two-input, two-output interferometer. Suppose we put three wells in a line 
as shown in Fig. 7a. Particles in the outer wells A and C are connected by tunneling to the center well 
B. The tunneling Hamiltonian, without interactions, is 

H/h= X{a^b +ab^+ b^c + bc^) (13) 

If we solve the equations for a{t), b{t), and c{t) we find that the solutions repeat the initial conditions after 
a time = 7 r/(-\/ 2 A) but at half that, tn = ^T/{2^/2X), the equations for the operators are equivalent to 
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Figure 6 : Effects of interactions 7 on the double-well HOM results with N = 8. (Cf. Fig. 3(Right)). 
(Left) Probabilities of particle occupation of the well A for initial conditions Na = Nb = 4, and 
interaction parameter 7 = 0.3. (Center) Similarly with 7 = 0.5. (Right) Similarly for 7 = 1 . 


detector equations for the interferometer shown in Fig. 7b. Initially we take a single particle in each of 
the outer wells. The initial configuration repeats in tu- If a Fock state for the occupation of the three 



Figure 7: (a) Three wells in a line, (b) The interferometer equivalent to the three-well line of (a). The 
phase shift needed is (/> = tt. 

wells is |mi, m 2 , TO 3 ), then the wave function at time tn becomes 

mn)) = |200) - ^ 1020) - ^ 1002) -f i |101) (14) 

Thus while the original state |I, 0,1) is present, the only other states are those with two particles in each 
well, in a kind of HOM effect. 

We can place four wells on a square as shown in Fig. 8 (a) with particles initially only in wells A and 
C. Tunneling is between nearest-neighbor wells with the initial state being {A,B,C,D} = {1010}. Possible 
states are (2000, 0200, 0020, 0002, 1100, 1010, 1001, 0110, 0101, 0011}. Now the repeat time for the 
initial wave function is 7 r/ 2 A. At half that, tn = 7 r/ 4 A, the system is equivalent to the interferometer 
shown. With a single particle in each of sites A and C, the wave function becomes 

^ (|2000) -f 10200) -f |0020) -f |0002)) -f i (|1010) - |0101)) (15) 

v8 ^ 

So in this case the interference removes occupation from these: {1100}, {1001}, {0110}, and {0011}, i.e., 
states with nearest neighbors both occupied. 
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Figure 8: (a) Four wells in a square with tunneling only along the sides. When all four tunneling rates 
have the same value, A, the system, at t = 7r/4A, is equivalent to interferometer shown in (b) with phase 
shifts C = 0 = 0. 


However we can do a more interesting experiment with this form of interferometer, namely look at 
a violation of the Bell theorem. An interferometer method of testing the Bell theorem with arbitrary 
numbers of bosons was proposed in Ref. [32], with Alice and Bob varying phase angles C and 0 in a device 
like that in Fig. 8(b). Here we consider the four-well system with just a single particle in each of sites A 
and C. However, in order to vary the phase angles in the four-well device, we have to allow for variable 
tunneling rates. Fig. 9 shows how we proceed and defines the variable tunneling rates Ai and A 2 . After a 



Figure 9: Four wells with Alice’s having counters A and D ; Bob has B and C. The tunneling rate between 
A and D is A as is that between B and C. The tunneling rate between A and B is Ai and that between 
C and D is A 2 . 

suitable time, t = tt/AX we measure the parity product {PaPp) with particles in A or C counting positive 
and those in B or D counting negative so Alice’s parity is Pa = (—1)"® and Bob’s is Pp = (—1)"^. We 
adjust Ai and A 2 to form situations P^ and P^ and maximize the usual quantity in the Bell-Clauser- 

Horne-Shimony-Holt (BCHSH) form [33] of the Bell theorem: Q = {PaPp) + {PaPp)pl^PaP'j^ — {PaP'^- 
Changing the tunneling rate is equivalent to a change of phase in an arm of the interferometer. The 
tunneling rate Ai is Alice’s setting and A 2 is Bob’s. Denote a parity average by ^). We maximized 

E{\, 1) -I- E{\ -I- 1) -I- E{1 ,1 — — E{1 + C) 1 ~ C) found the maximum was 2.815 at ^ = 2.74. This 






result is very close to the maximally possible violation of 2-\/2 [34]. 


5 Conclusion 

We have shown how a simple set of potential-well traps can be used to demonstrate the HOM effect and 
to build interferometers where the beam splitter is provided by tunneling between wells. This suggested 
method emphasizes the possibility of doing generalized HOM experiments with multiple particles by 
using condensates. The actual implementation of such interference effects involves some requirements 
that could be difficult: initially counting particle numbers in the wells when used as sources, turning on 
tunneling for a specified time, and then counting the particle numbers in the wells acting as detectors. 
Ref. [11] has shown hat these difficulties can be overcome in the N = 2 case. When interactions among 
the particles are considered we find the HOM effects change considerably, showing the effects of possible 
“fermionization”, cooperative pair tunneling, and self-trapping that have been seen in previous double-well 
analyses. We find that rearranging the wells in various configurations can allow their use as a wide range 
of interferometers. 

After completing our work we learned of related calculations by B. Gertjerenken and P. G. Kevrekidis 
[35]. We thank them for sharing a preprint with us. We thank Dr. Asaad Sakhel for a critical reading of 
the manuscript. 


Appendix. Self-trapping 


We referred to self-trapping as a possible cause of the results shown in Fig. 6. We would like to examine 
that claim a bit in this appendix. Ref. [28] discusses the dynamics of a Bose-Einstein condensate in a 
double-well potential and the “self trapping” effect, which occurs when the interactions between particles 
change the characteristics of the tunneling between the wells. Within a mean-field approximation, free 
tunneling between wells is reduced for interactions above a certain critical value, causing an initial im¬ 
balance between numbers of particles in the wells to be locked in with only smaller amplitude oscillations 
away from that. In this approximation, all particles occupy the same state given by a Gross-Pitaevskii 
equation 

■i/'(r) = fci'Ui(r) -I- k 2 Ui{r) (16) 

where Ui{r) represents the groundstate wave function in well i = 1,2. The mean-field equations are 
expressions for the amplitudes ki,. When a more exact analysis is used (a full N-body wave function 
depending on more than two parameters, but still a two-mode approximation), quantum fluctuations 
alter the mean-field results by time modulations and system revivals. 

If we derive equations of motion for the annhilation operators in the Hamiltonian of Eq. (2), we find 


.da 


i- 

i 


dt 

dt 


—b + "fa^aa 
—a + jb^bb 


(17) 


where we have measured energy from Eq = 0 and time in units of h/X as before. If we replace a^a by 
NA{t) and b^b by Nsit) these are just the mean-field equations for the ki found in Ref. [28], which, if 
one has initial conditions AfA(O) = N, have the solution [28] 


NAit) = y [1 + cn(2t|Af^7^/16)] 


(18) 
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where cn(u|z;) is a Jacobi elliptic function. An equivalent set of equations suitable for numerical solution 
for arbitrary initial conditions is given in Ref. [27]. We can compare our own exact results of Eqs. (7) 
with these mean-field results. 

Note that in Sec. 3 our initial condition was always Na{0) = Nb{0). In this case one finds in both 
exact and mean-field cases that NA{t) = N/2 because state {n,N — n} is always paired in a superposition 
with state {N — n,n}. (Note the symmetry of the probabilities in Fig. 6.) However, self-trapping shows 
up dramatically when all the particles are initially in one well, Afyi(O) = N. The critical interaction for 
fixed N is 



(19) 


In Fig. 10 we consider a small value N = 8, although our equations are valid for any N value. The critical 
value then is jc = 0.5. Our plots show the time dependence of iVyi(t). For 7 = 0 the mean-field and the 
exact value coincide with a sinusoidal oscillation with full amplitude 0 < AOi(t) < N. At 7 = 0.3 < 7c the 
mean-field continues to show full oscillation but an additional fluctuation is introduced in the exact result. 
At critical value 7 = 0.5 the mean-field expression reaches self-trapping with NA{t) = N/2-, the exact 
result fluctuates about this. (This result favoring equal populations in both wells might be relevant in our 
own case of Fig. 6 .) In the case of 7 > 7c the mean-field result maintains a value of A^yi(t) > N/2 with 
oscillation amplitude decreasing with increasing 7; the exact quantum value seems to have an additional 
longer period oscillation frequency superposed, similar to the two-frequency oscillations discussed in Ref. 
[24]. If Na{0) < N the amplitude of the oscillation below criticality decreases correspondingly around 
N/2. 

None of these results exactly mirrors what we found in Fig. 6 , where the configuration {4,4} became 
favored. In that case, the initial condition was Na{0) = 4, for which the symmetry of the probabilites 
leads to NaII) = 4, i.e., constant in value, both exactly and in mean-field approximation. However, if 
we look at the details of the probabilities Pn{t) = |cn(t)| for individual configurations (n, N — n}, then 
we find that the initial configuration itself can become “trapped” at a large value, as seen in Fig. 11 
for a large interaction parameter. This effect is analogous to the “fermionization” we saw for just two 
particles in Fig. 5. This effect is due simply to fact that the initial condition enters with approximately 
unit amplitude and that any other states are mixed in only in higher order in I/ 7 . 

This is an effect that cannot be seen in the mean-held calculations; we plan to return to this subject 
in a future publication. 
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0 5 10 IS 20 25 ‘ 0 5 10 IS 20 25 ‘ 

Figure 10: Time dependence of NA{t) in a double-well potential with N = 8 with initial condition 
Na{0) = N. Solid line, exact result; dashed line, mean-field result of Eq. (18). (Upper Left) No 
interaction 7 = 0 . The exact and mean-field results coincide. (Upper Right) 7 = 0.3, below the self¬ 
trapping critical interaction value. (Lower Left) 7 = 0.5 at the critical value. (Lower Right) 7 = 1.0, 
above the critical value. 


kn(t)P 



Figure 11: Time dependence of configurations probabilities |c„(t)|^ in a double-well potential with N = 8 
and initial condition Na{0) = 4 with inteaction 7 = 10. The solid line gives |c 4 (t)| for initial configuration 
{4,4}, while the dashed line is for |c 3 (t)| = |c 5 (t)| for states {3,5} and {5,3} respectively. All other 
states have very small probabilities. The {4,4} state is trapped at a relatively large values. 
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